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ABSTRACT 

For any system whose input-output relationship may be 
described by a ratio of polynomials in S (the complex number Vtyw), 
one of the measures of the system performance is its transient 
response to astep input. For systems of order greater than second 
the determination of the parameters of this transient response 
has always been a difficult process, In this thesis a new 
procedure is developed which allows the rapid and accurate 
determination of the transient parameters of any system up to 
and including third order with one zero, The procedure consists 
of developing the response equation in terms of angles at the 
Singularities rather than the magnitudes. The resulting equation 
is utilized to develop expressions for times of extremals, 
magnitudes of extremals, and a form of rise time. A simple but 
accurate graphical method is demonstrated for solving the result- 
ing transcendental equations. It is shown that the expression for 
times of extremals is a function of only two variables and may 
therefore be displayed on a two-dimensional plot. The times 
determined from this plot are utilized in computing the magnitudes 
of extremals. Various curves are presented to provide a means for 


general analysis. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


Definition 
The complex frequency (S=V4y), 


The real part of a complex number; 
neper frequency. 


The imaginary part of a complex 
number 5 radian frequency. 


The S-plane coordinates of a pair of 
complex poles. 


The natural frequency (l= V7 Y* +a), 


The S=-plane coordinate of a real zero. 
The S-plane coordinate of a real pole. 


The angle measured clock wise from 
the negative real axis to the vector 
from the origin to the complex pole 


Va td w,/. 


The angle measured counter clock wise 
from the real axis to the vector from 
a real pole to the complex pole Td tg ats. 


The angle measured counter clock wise 
from the real axis to the vector from 

a real zero to the complex pole Val tg ats 
Real time, 


Nondimensionalized time (7 =4/), 


The time to the nb point of zero slope 
of the transient response. 


The time to the n& point at which the 
transient response is equal to the 
final value. 


The time response, 


Normalized magnitude of the transient 
response at time ton (Mpdn = Xion) Jbinol Valve). 


lion Jintroduction. 

In recent years the study of systems which have an input- 
output relationship (or a relationship between any two variables 
of a system) which is described by a differential equation has 
increased tremendously. Such relationships are found in the 
equations of motion of objects, in the equations for chemical 
reactions, in thermodynamic relations, in the equations for heat 
and fluid flow, in electrical network equations, in automatic 
control and servomechanism equations, and even in the equations 
governing the fields of industry and management. These relation- 
ships are written in either of two ways: 

(a) Given that the variables of interest are X(t) and Y(t) 

which are expressed as functions of some third variable 


W(t) and its derivatives: 


nel ne 
Xd) = an OW YI + Tyo! A wh) SE a a ee tao Ww 


m mel 
Yu)= Pre ML + bos? yr oe . - e  ee 
(b) Given that the variables of interest are X(t) and Y(t) 


which are expressed implicitly in the equation: 
n 
2, 7 Myr + -- i — thx tan? Fars ----- t4ey 


In either casé the complex frequency (S = vaya) may be introduced 


and the relationship expressed as¢ 





(1) Xs) — Ans"+ nei S "4 --- - 7 7 7 a 
Yo bins + Cen + ae a Sie nee + bo 


The usual case of interest is to allow Y(t) to be some input 
(either an initial condition or a specific function) and to 


investigate the response of X(t). A specific case which is often 
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mentioned in design specifications is the response of X(t) when 
Y(t) is a unit step input (i.e., a unit value input applied 
instantaneously at time zero). Depending on the order of the 
relationship , the value of Min the equation (1), the determinat- 
ion of the response may be straightforward or quite tedious and 
time consuming. For the pure second-order system (m=2 and n=0) 
the response has been extensively investigated, and general curves 
and equations are available in any elementary text on servo- 
mechanisms. For relationships with m2 and/or n?0 the information 
is not nearly so plentiful. The approach taken by most authors 
is to present the equation in basic form, comment briefly on the 
effects of the higher order terms, and turn to methods of making 
these higher order terms negligible in order that the response may 
be analyzed in terms of the more familiar second-order relation- 
Ships. Unfortunately, the higher order terms cannot always be 
made negligible, and the engineer's only recourse has been to work 
out each system in detail. No method has been available for a 
general analysis of the effects of the higher order terms. In this 
thesis a procedure is developed which enables the analyst or designer 
to solve in general an equation of one higher order in both 
numerator and denominator (i.e., order three in denominator and 
order one in numerator.) 

Perhaps the simplest way to gain an appreciation of the problem 
is to examine the relationship on the S-plane. Equation (1) may 


be rewritten as follows: 


XiMer _(srelsepy - --- - Gr %) 
Vis) (SAA) (S92) - ~-—~ (Stfm) 


Where the Zi and Pi may have either or both a real part and an 


J (a) 


imaginary part. 


An example of such a system appears in Figure l. 





Figure ls: An example of S-plane diagramming 


The system depicted would be written as: 
ae (31 8NSi 22) 
Sire $(3+Pi) (54 Ps) (S47, Je )GIT, Joye Tat quid V2 <Jwz) 
where the S term in the denominator represents Pin, the step input. 





The usual nomenclature assigns the term poles to the P 
values (represented by D ) and the term zeroes to the Z values 
(represented by 0) with the collective name of singularities assign- 
ed to all roots of both the numerator and denominator polynomials. 
These singularities may occur anywhere in the S-plane and in any 
combination, but for stable physical systems the poles always occur 
in the left half plane (T<o), and complex singularities always 
occur in conjugate pairs. 

The second-order approximation approach is to consider the 
"dominant" pair of complex poles (9; $y, in the @xample) as the 


determining factor and force all other singularities far out in 





the left half plane so that they may be considered negligible. 

A generalized method of analysis for systems with one pair 
of complex poles and one real singularity has been developed 
by Patton and Abbott [1]. This method consists of using a 
scale factor to transform the real singularity to the V=-d/ 
‘position on the real axis and then determine various response 
parameters by entering a set of generalized curves with the 
location of the scaled complex poles in relation to the fixed 
Singularity. A study of this method was the starting point in 
the search for a method of generalizing the equation for a system 


with a pair of complex poles, one real pole, and one real zero. 





2. Attempts at Generalizing the Equation for a Third Order 

System with One Zero. 

A clarification should be made of what is meant by general- 
izing the equation. If the stated problem contains more than 
two variables, it is desirable to combine the variables by means 
of some linear transformation so that a maximum of three trans- 
formed quantities remain as variables and all other variables 
are "buried" in the transformation, 

This point is illustrated by means of two examples: 

(1) The pure second-order response to a unit step input 

is expressed as: 

Xes) * rae aa 
Assume that a time (time of peak, rise time, settling 
time, etc.) is the transient parameter of interest. 
There are then three variables of interest: Vel, ab, € ¢ 
(k appears only as a multiplier). It is possible to 
display the information desired on a two-dimensional 
graph by displaying Y on the abscissa, t on the 
ordinate, and “/ as a parameter for a family of curves; 
however ,the problem may be simplified by noting that 


the time solution may be written ass 
Se w 
Xay= wf [1- We the e Sin (wd +) | 
Where T= itn fe “k,) 


Now let T= yt , and U= Tah : 
Xue ee U PE ee ae (T+r)] 


Ce Airc Laps ( Ka ) 
This information could be displayed on a two-dimensional 
graph with U as abcissa and T as ordinate. The necessity 
for the family of curves has been precluded by "burying" 
the variable @//in the transformation t= uyt ; 
(2) As a second example, consider the response of the pure third- 
order system to a unit step input: 
Xs) = S57 (s+ 2WSF On) 
There are now four variables: G/ w/, 4 fe 2. Obviously, four 


variables cannot be displayed on a two dimensional plane. 


Writing os response en the © mee ; domain: of 
10% - =f ?¢ uJ? (7 y . 
XT =o #,[1- {We Pgtd® PE — sin (ugh + 7-07?) 
N* iin? (Eton ta) wel hij? + (P- Tol) 


Wl = pele Um) 5 Op = predan ( “pf. i) 
Now let Taubdd U= Yy/ and V= Aur/: 


Lt Cee ee 7 
X (Yd) = a a7? S27 Z+Wv- w)* te ri (TH+ aye 


i= fred (44) > Ge= Acetan a) 


This information may be displayed on a two-dimensional graph 
by plotting U on the abcissa, T on the ordinate and V as a 
parameter for a family curves. Again, the fourth variable has 
been "buried" in the transformation T=uyyf ; 

It may be noted that the transformation used in the above examples 


was not the only possibility. In either examples, another 


possibility would have been T=Vt; and in the case of the third 
order system, another possibility would have been T=pt. The 
latter is the transformation used by Patton and Abbott c7] who 
showed that this transformation is equivalent to placing the pole 
at minus one on a W-plane defined by W=s/p. These authors then 
displayed the response parameters of rise time, time of peak 
overshoot, magnitude of peak overshoot, and settling time by 
providing curves on the W-plane representing constant values of 
these parameters. The parameter values for a particular system 
may be determined by locating the conjugate roots on the W-plane 
and reading the value of the curve passing through these points, 

The response of a third-order system with one zero to a unit 
step input is expressed as: 

Ke K0$7-2) 
X (8) ~ S(StP)(S*42TAS + Wp) 

The number of variables has now been increased to five: wy, U/,%4P, fz, 
This number can be reduced to four by means of a transformation 
Similar to the type just described, but a maximum of three may be 
displayed on a two-dimensional curve. In order to reduce the 
problem a transformation must be found which "buries" two of the 
variables. One approach is made by noting that the pure third 
order system was generalized by making a transformation which 
placed the real pole at the minus one position on a W-plane,. 
It appears that if a transformation could be found which would 


fix the two real singularities at two predetermined points on some 


W-plane where W is a linear function of 5 then the only variables 
would be the angle and distance of the conjugate poles from the 
origin of this W plane. In order for this transformation to be 
useful there must exist a unique relationship between the 
fictitious time associated with the W-plane and real time associated 
with the S-plane, and likewise such a relationship must exist 
between fictitious magnitude and real magnitude. 
Although many minor attempts were made to determine a trans- 
formation of this type, only two seemed to have any merit: 
(1) The bilinear transformation. 
It is well known that by means of the bilinear transfor- 
mation any three points on the S-plane can be mapped 
into any three points on the W-plane, If it is desired 
to map the points P, Z, and origin into the points -l, 


-2, and the origin, this may be achieved by the 


transformation: 
(Wt1)(-2) _ G-AzZ 
W [-2- -1)] ae -) 


which reduces tos 


Vee es = 6> Jee 
male) S igen W (2P-Z)42(P- Z) 


(2) Transformation by scaling and translation. 
Another means of placing the pole at the minus one 
point and the zero at the minus two point is by 
compressing (or expanding) the scale until the distance 


from the pole to the minus one point is equal to the 
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distance from the zero to the minus two point and 
then translating through this distance. This may be 


accomplished by the transformations: 
we (eR) [st (2-2°)] 
5s = (p-2) [w+ GE2)) 
Although either of these transformations can be shown to 
be suitable for certain specific cases no unique relationships 
could be determined which would apply for the general case. 
Nofurther attempts were made at determining a trans- 
formation, Although no proof is given, it will be shown in the 
next section that it is unlikely that a suitable transformation 
exists. 
One of the major difficulties experienced by Patton and Abbott 
even in view of the generalizing transformtion, was the solution 


of a transcendental equation of the form: 
- aw 


ke SS 7 (aly PED) 
In computing values for the curves in their thesis they were forced 
to use ‘iterative methods for the solution, While studying their 
work a graphical method was discovered for solving this form 
of equation, and this discovery led to the idea for a means of 


generalizing the third order system, 
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3. Development of the General Equation. 
The response of the third order system with one zero is 


expressed in the 5 domain as: 





es k (S72 L Sea ee 
X¢s)* S(SIP)(S?¢2TyStitn) — SCS 1P) (89 Tol 4g wel USF L/- Ju) 
and is expressed in the time domain as: pt 





eZ ee ee 
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oppaall ke OE § Leer o) 2 BLA Cos uy 
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The time domain equation is derived by simple residue methods and 


algebraic manipulation, The various values are depicted in 


h yw 


Figure 2, 





Sie 


Figure 2: S-plane plot of noninnnnnnneesaine 
S(54P)(S 7424S 49474 w/?) 





The expression for X(t) may be simplified by rewriting it in 


LO 











The following relationships 


terms of the angles @p, @z, and v. 


are written by inspection from Fig. 2: 
Vs @ Cy” st won 

PD“ = (P-TA “+ u}/* 

a s 2 a 

Bp = (2-wh + a/ 
SN = Wl Cos r= hon 
Sin@ps wel/ Pp; casas P- Wel) DB 
ud/Zp 3, CosOg = 2- WS, 


Sindg = 


By trigonometric identities: 
Sin @2-Ha) = He Ae) 

2 

COS (@2-Gp) = E-WAE Ne) tue 


20 PD 
a 
Sin (ult + G2 - spk 


From these relationships: J = 
Wp2d s; " can PZ) € 
Mee Sin (G2 -Gaje FD sinv 





Xa) = ke s1+ 
@) Wy? P 2 Sn PO 
By further use of the relationships above: 
an 31 Of 
Wn Sindg sine (Cet Og F Ce tr) 





Wp ZD ~ wr WK Op bane 
& PD wh (604 Oz +E2 10) se, 
Sia _ OW Ral ew ae 
Siri (827 a) 


CoSOB SING 45GB SH 
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Sot 


na J 
PZ0_ _ wel (cotop Co /r) WL oe , 
wel (colbz +co/2) ui fe 9/2 Sra 


ZPD 
= S10 (opt) 
~ 6,7 (62 FOO 
-[eo/ep tca wae 


I Ef Si n(6e-%2) oe 


Therefore: 
SAT S/ACCe +7) y 
) 9) ey 
| 482 -G mf 


Eg'n. Te Xt) = £é, jf r 
_ Sintepin § S17 (ul 


SIF Sen (O29) 
Eq'ne I has a number of important properties which are noted here: 


aft 


(1) For stable systems the multiplying constant IRS is the 
final value determined bys 
Final Value = see JS Xos) 
SO 
(2) If the transformation T=4)/ is made the equation is 
completely non-dimensionalized; i.e., the equation is 
a function of angles only. 
(3) The equation may be normalized to unity by multiplying 
by ae work 3 ieee the final value has been 


normalized to unity. Using these properties: 


’ ») 
: t os SING 3/4 - — te +Cor@ 7 
Eg'n, dae Xr (Ga) = y $ 5 3 a SF ) (co / 


in (p10? ~coIOT 
> sin elr . Sin (1402-07 +4) 


SF Sin (arr 


Eq'n Ia is valid for @z#180- % and 6p < 180"= % ¢ i.e. for the 
zero not at the origin and the pole in the left half plane (it 
is assumed that 7 92). When Je) 2180-9 the system is unstable 
and the final value theorem does not hold. When @z=180™ 2” the 
zero is at the origin and the final value is zero, These special 
cases are treated in detail in later sections, 

Eq'n Ia is valid for any order system up to and including third 
order with one zero, Any combination may be achieved by allow- 
ing either or both of the angles Op and 62 to go to zero which is 
equivalent to locating the singularity at infinity, 
This equation may now be manipulated to determine various parameters 


of the time response, 


aw 





4. Development of the Equations for the Parameters of the Time 

Response. 

For a given system the slope of the response curve at time 
zero can be quickly computed along with the final value of the 
response. Since it is also known that any oscillatory motion 
is of the damped sinusoid form it is necessary to know only three 
other parameters to be able to sketch quite accurately the 
complete transient response. These parameters are: 

(1) tc-time at which the response crosses the final value. 

(2) tp-time of extremals (points where slope of time response 

is zero.) 


(3) Mpt-normalized magnitude of the response at the extremals. 


[Mp4 = Xn (4e) 


(a) Development of equation for the time of extremals (tp): 
In order to determine this value the expression for 


X(t) is differentiated and equated to zero: 


Ut [Xen] * EB, } (Ses feet) ep) 
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equating to zero: 


13 





There is no "closed" method for solving this equation for tp, 
but a simple graphical method exists and will be explained 
presently. Note that the times of the extremals are independent 
of the angle 0. 

(b) Developments of the equation for normalized magnitude of 

extremals (Mpt): To determine this value the equation for 
X(t) is entered with the results of part (a): 


Mpt = Xunfbee =77 a Sin OP Sin (O2-G Cee 


Sin TT Sin CO271T7) 


~o% 
~ LANBID). © See (w}pe+82 -92 747) 


SF $'7(C270) 


Substituting from Eq'n ree the quantity in braces: 
Agn.m: Mpt = 1 p SilGe7. 2- Worl sino? Sin(umlp eG) Sialalprerepr)f 


Sn TSN lO2ZIT) Sin (p77) 


Values for solution of this equation are determined graphically 
at the same time Eq'n II is solved, 
(c) Development of the equation for the times when the time 
response crosses the final value (tc): 
This time occurs when the last two terms of Eq'n I are 


equated to zero: 


7 She 1h (029) ~The 
Sin & Sin G2 Lp) = oe Se S77 (ol fe 482 -O7 +7) 
SAV 5+ (6277) S19 Sinl(C277) 


op -tly, 
Egan TL 22S We? 5. (ute 102-Cp tT) 


S10(6P17) 
Again tc may be determined graphically on the same plot used to 
solve Eq'ns II and III. 


The notation of the equations is simplified by making the 

: JS wt 
transformations as before: We 
kr, = “Wns 


1 





The transformed equations then becomes 


- Pp, 
kgr. La » Sin (2-Gp) 2 ee Sin (7p + @2-8) 
‘ e -7 Or) 
Egn. la. Xn (B= Mp -}1 7+ S17 (Ett). 2 Ze "[ Sin 8 Sin (19762-9) 


SPT S) (9292) 5,2 (opto) 


—- Sin (Ip + 2-P olf 
a =e Sin (62-&P) = % fy 3b 
$ se Sin (Te 402-9217) 
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S. Graphical Solution of the Equations. 
Equations Ila and [Va may be solved graphically to an accuracy 
of approximately one percent by following the step by step 
procedure outlined below. Equation IIIa may be solved directly after 
the solution of Equation IIa for Tp. 
(a) Graphical solution is performed on semi-log paper with 
values of T plotted on the linear axis and values of the 
functions of T plotted on the logarithmic axis. 
(b) Compute or measure the angles @p, 9z, and? on the S-plane 
diagram of the problem to be solved. 


(c) Compute the following values: 
Sin (@z-8p) 
Sin ep/Sin (@p+7) 


(d) If Sin(@z-@p) > O Construct the function Sin T. 
If Sin(@z-@p)< O Construct the function Sin(-T). 
(Note: This function is easily constructed by plotting 
the values of Sin T for T= 10°, 20°, 30°, 5°, 
60°, 90° and appropriate multiples. For T < 10° 
the function is approximately a straight line.) 
(e) Construct the function Sin(@z-@p) e “ “P 
(Note: On the semi-log paper this function is a straight 
line with intercept Sin(6z-@p) and slope of - as Fs +) 
(f) Locate the point on the straight line where the directed 


horizontal distance to the sine function is equal to 


(6z-Op). The value of T at this point is Tp. 
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in Op $10 (9g- 9p) “Tha o 
(g) Construct the function $10 Ge Soap nee Se ~~ 


Sim (9277) 


(Notes If iq Se 72) is of the opposite sign to 


Sin(@z-@p) it will be necessary to make a second 
plot on semi-log paper following the rules laid 
down in (d) above.) 
(h) Locate the point on the second straight line where the 
directed horizontal distance to the sine function is 
equal to (@z-Op+ 7). The value of T at this point is Tc. 
(i) Compute the value of Mpt by use of Equation IIa and the 
value of Tp. 


(j) Use the inverse transformations: 


f= T el 
Xp)? Ee Mps = pie Met 


to obtain true parameter values. 





(k) Repeat steps (f), (h), and (i) for as many values as 
necessary to obtain a complete description of the transient 
response. 

The crucial step in the graphical solution just outlined is the 
determination of the sign placed on the values of the logarithmic 
axis, Since the construction is performed _on semi-log paper it is 
impossible to show both positive and negative values on the same 
plots; consequently, only alternate half-cycles of ot may be 
nroe eed. If 09z2< Op the vameror sin(@z-@8p) is negative and the 
logarithmic axis must be labelled with negative values; therefore, 


the negative portion of SinT must be plotted. Likewise, if @z > 6p 
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the value of Sin(@z-6p) is positive and positive portions 
of SinT must be plotted. Similar rules apply to the solution 
of Equation IVa. 
The procedure is illustrated by means of the following 
example: 


Given a system for whichs: 


On = 60° Jl = 1 
O27 =40% kz a7 
y= ey ° Wn *P 


The computation proceeds as follows (see Fig. 3): 

(1) Construct the function SinT on semi-log paper. 
(Note: 6z-@p<O; therefore, plot the negative 

half-cycles). 

(2) Construct the line Sin(@z-@p) 27 Hen SF 
(Notes Qz-Op<o; therefore, Sin(@z-@p)< 0) 

(3) The equation to be solved is: 

—,_ 643 2” Veith Sin (T-40%) 

Locate the points where the horizontal distance 
between the straight line and the sine function plot 
is )0° (straight line point to the right of the sine 


function point). The T axis value of the points on 


the straight line are the values Tp: 
Tp = 4 1 radiens 
7Tp2 = Lowe radhans 
TPs = /0.3 radians 
Note that these values do not depend on the value of 


0 but only on the relative location of the complex 


poles and the real singularities, 
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Fig.3: Computational Diagram for the 
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(4) 


(5) 


(6) 


Construct the line; 


Sin (Get) oo Veor 


Sine SinlO@2r7r) 
corn! 7/4. Ss 


= .335 C 
Read off the values of this line at the times Tpn: 
(Call the values Cn.) 
C, = 213 
C,, = 165 
C3 = 113 
(Notes These values are positive since the logarithmic 
axis was assumed positive for construction of 
the line.) 
Read the values from the sine function at the 
correspondence points found in step (3): (Call 
these values Bn.) 
Bj=- son 
E5=s at? 
By=- 0103 
These are the values of Sin(Tpnt6z-@p). 
Add the distance (0z-@p+%)=l° to each value of Tpn 
and read the corresponding sine value: (call the 


value An.) 


A,=- 1.0 
An== 1.0 


A3= Sinh33°=sin73°=.955 


These are the values Sin(Tpn+@z-9ptT). 
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(7) Compute the normalized magnitudes of the extremals: 


- 517 (0e4 ym) “Pont Sine a - 
Mpdn * pr = S08? Sin (Tent Gg - 9p) Sin (Tort 8g 67T) 


HAY SinlG¢+7) Sin (opt) 


= 7 ¢ Cn [3.58 Ba-An J 


Mpt3 =] ,0715 ; 
(8) Construct the line Js 
= n@ 


Sin (opt) 
a T/S3ee 
= «Zee 7 
(9) The equation to be solved is: 
aitene oC 
~£.3€ . = Sin (TH44%) 


Proceed as in step (3) to find these points of correspondence, 
The value of T at these points corresponds to the times when 


the response crosses the final value (Tc): 


Tc, = 9,1 radians 
Tezg.=11.6 radians 


This completes the values necessary. 


c | a Se 











Figure h: Step response for O;7= 90°; Gp= 50 * Ts oo a/= 4 
In order to demonstrate the accuracy of the method, this 


system was simulated on the CDC160) digital computer using a 
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method of solution completely independent of the results obtained 
here, The computer response is shown in Figure 5 (Note:@d on the 
computer Simulation was 2.0). For purposes of comparison, the 


computed values and the computer response values are tabulated: 


Parameter Computed Computer Response Difference 


Tp, leat 4.3 6.8% 
ips 6.82 6.9 2 5% 
TDs 1033 10.3 a 96% 
Mpt, ® 9766 @ 968 1.2% 
Mptz foe 728 055% 
Mpts 1.0715 1.076 51% 
tc, yet oe 1.1% 
TCz 1436 Pie 1.75% 


(Notes Computer response values were interpolated from 
the computer readout tabulation given in 0.1 


second steps) 
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rig.5: Computer Response for the 
: G <? 0 
System: 9271,0;6p-80; 7-44. 
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6. Adaptation of the Equations to Other Singularity Combinations, 

It has been stated that the equations developed are general 
in that they are correct for pure second-order systems, second- 
order systems with one zero, pure third order systems, and third 
order systems with one zero. This property is developed in 
more detail in this and later sections. Also, although an 
important value of the work just described is the method of solv- 
ing the transcendental equations to obtain various parameters of 
the time response, an equally important value is that the means 
have been provided to develop curves and relationships which 
allow the analyst and/or designer to observe the effects of the 
real singularities at a glance. This point is also developed in 
detail in later sections, 

To make equations I through IV applicable to the particular 
system of concern, it is necessary to modify the multiplying 
constant > £2 Japp « it is more correct to define the 
multiplying constant as the "Final Value" of the system. The 
final value theorem for stable systems states: 

Je Say) = im S Fes) 
} — 0d SPO 


The final values of the various systems are tabulated in Table 


1 along with the values of the angles 92 and Op. 


stem Order'| Zere Final Value | @2 | 6 


Two None K/wry® Zero Zero 
Two One k2/wn* Variable Zero 
Three None k wr? ? Zero Variable 
Three One k2/v,?P Variable Variable 


Table 1: Modifying Relationships for Equations I through IV. 
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As an example equations I through IV (modified) are listed 
for order two with no zeros (pure second order): 


« (coda) ual? 
Eqn. Z (mod): X41) = # 5 -~,¢ S10 (ul) +r)f 


0 = Sin wtp = Wdippa= NT 


= (cody) nT 


Egn. DT (moe): 
Cos niy 


Mp 4,.= i-@ 


Egin ID (red): 
os Son (let T) > Why Fe nT- 


Eg'n If (med): 


These equations are well known for the second order system 


although the form may be different. 
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ie Analysis of the second order system with one zero. 


Following Table 1 of the ~— section: 


© (code ups — 
Eg\n. ,LTimeod ); Xa) £2 241-545 . Sin (ul 4 Be 120)5 


£gn Zim? 02 Sin (welt 403) => wilipn= yy -O2 
\ e e elriudipn 
£gn. Kun); Ma 4 2 pi- 8 air isin (leila per 


=$i- GaomeeT = | sin 5 


$'9(G3 17) 


Eg'n. TW invl)s OF SAV tgs 10217) = dplen = NIT - O32 - 
which are the governing equations of the second order system with 
one zero. It may be noted that each of the time response parameters 
may be expressed as the pure second order parameter with a 
correction factor. In the case of the times of extremals and 
times when response crosses the final value, the correction is a 
linear one; i.e., the time is decreased by an amount O4p,/, The 
correction factor for the magnitude of extremals is a multiplicative 
nonlinear factor. The magnitude of extremals may be written: 

X(tp) = Final Value (1 + % overshoot) 


-(eolnT 
ands % overshoot (second order) = = COS(n7)¢ 


[Sing 


~ (ode) [nr-bal 
SintP err) c 


and: % overshoot (second order with one zero) = - COS/NT [—— 


- (eodr) b3 


= | Zoversheod trecenclorckees] $03 ae 
the correction factor =| 3147 6 bzcot¥ 
=nl@217) 
Figure 6 is a plot of Mpt of a second order system with one 
zero versus COS@for various values of @z. The curve for 6z=0 is 
a plot of the percent overshoot of a pure second order system, 


The system fails when 02=77-7% 3; i.e. when the zero is at the 


origin of the S-plane. At this point the final value is zero but 
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computed percent overshoot becomes infinite resulting in an 
undefined value. It is necessary to return to magnitudes to 


interpret this special case: ; 
- Co x (nt-02) 
Xe) = ae LE Cosnp 2 e 


Wy, + Sin (O27) 
rz J _codr(nt-62) 
- —_—_ = e Ww 


a ZO lyn 
fe . kcosn7 = 22 a7 cdtr (nt-32) 


Uy, 4 2 


be Hing G2= 1-7 5 2D ty 5 ZO 


\\ 


x = Gof to-1) 7 eae 
(+p) jee ae ae C 


for the curves of Figure 6,n=1 . 


, = kee 
a) Xp]... a Cn C 


The equations obtained here are quite simple to use and no 


further analysis of the system will be undertaken in this thesis. 
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Second Order System with One Zero 
(Final value normalized to unity). 
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Fig.6: Magnitude of First Extremal of a 
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8, Analysis of the pure third order system. 


Following Table 1 of Section 6: 
- [cote eed schyail 


(Sin G/ 
Eg*n. Z (med)? Xone U3) as te (222 = S26 21 5 in (alll 40-6 


- wtp) for 9p - 


Egn. I (mol)? - Sinape iat 


- j on 
Eg’n. (med), pt = 31 = oe call i‘ Sin (ube -bp) = Sin bade or-apye 
2 ae, 
Eg'n, Bt (we): -=27% ¢ a Shrt? = 5 (wy, *7-Gc) 
Sin (ott) 


It is obvious that no simple correction factor exists for the 
parameters as it did in the case of the second order system with 
one zero. The methods developed in section 5 have been utilized 
to develop curves for the parameters of the time response. Figure 7 
is a plot of Tp versus Op (Note that Tp is independent of the 
enews 0). The curve has been extended to provide times for the 
occurence of the first three extremals (points of zero slope). 
It is interesting to note that the time correction is practically 
linear to @p=70° and may be approximately expressed as Tp=n7+ 6p up 
to that point. It will be shown presently that no overshoot can 
exist for @p 790°. For purposes of comparison the times of extremals 
are shown for a second order system with one zero (dashed lines). 
Figure 8(a) is a plot of the normalized magnitude of the first 
extremal of the response as a function of the cosine of 0 for 
various values of 6p (final value=1). Figures 8(b) and 8(c) are 
the same type of plot, but for the normalized magnitudes of the 
second and third extremals, The curve entitled 8p=0 is the 


familiar second order curve. An important point should be noted here, 
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For a second order system the first overshoot is always the 
greatest; but this is not necessarily the case for a third order 
system, For @p>76°, and dependent on the value of ¥, the first 
"overshoot" may actually be less than the final value (in such 

a case the term "overshoot" is misleading); for Bp>5°, and 
dependent on the value of 7, the second "overshoot" may be 

greater than the first. These effects are readily apparant upon 
examination of Figures 8(a) and 8(c). An interesting point is the 
case of 8p>90° ; on the semi-log paper used for computational 


purposes the line defined by the following equations 
ye Sin OPEC - Tiron ep 


is a line of positive slope and can never intersect the curve 
defined by the equation: 

y = - Sin T 
therefore, in accordance with the methods developed in Section 5, 
no time of overshoot can occur. This implies that for @p>90°, 
regardless of the value of 0 , the time response is monotonically 
increasing and never exceeds the final value. Figure 9 is a plot 
of te for the third order system (the rise time, tc, is defined 
to be the first time when the response is exactly equal to the 
final value). For values of @p>76° it is possible, for certain 
values of 0, that the response does not cross the final value until 
after the first peak. It should be remarked here that the curves 
presented in Figures 7-9 are intended to serve as a means of general 
analysis and not necessarily as computational aids. The analyst 
or designer interested in exact values should perform the calculations 
as outlined in Section 5. 
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An interesting point in the study of the third order 
system is the value @p=76°. This value appears to be the optimum 
setting to minimize fluctuations around the final value regardless 
of the value of 7 ; the response always reaches the final value on 


first rise and the first undershoot is minimum for all J. 
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9. Analysis of the Third Order System with One Zero. 

The character of the response of the third order system 
with one zero is best analyzed by plotting the values of Tp as 
shown in section 5, this information may be depicted on a two- 
dimensional plot since Tp is dependent only on the values of 6z 
and 6p. This plot may be constructed rather quickly and easily with 
Equation II by fixing @p (which fixes the slope of the straight 
line plot on the semi-log paper) and determining Tp for various 
values of @z. Curves of this form are displayed in Figures 10(a), 
10(b) and 10(c) for Tp corresponding to the first three extremals 
of the response. If the values of 6p and 6z are known these 
curves may be entered and the value of Tp for each of the first 
three extremals may be determined. Once this value is known, the 
magnitudes of the peaks may be quickly computed from Equation III. 
Again, it is emphasized that the accuracy on the curves is only 
about * 5%. If accuracy is required, the computational methods 
of section 5 should be utilized, 

A great deal of knowledge of the behavior of the system 
response can be determined from a study of the curves of Figure 10, 
This information is summarized in part in Figure 11. The coordinate 
plane used in Figure 10 is depicted again in Figure 11, but the 
Tp curves are deleted and the plane is divided into regions as 


follows: 
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These regions are formed by the line ep=90°, which is the angle 


at which the real pole moves "inside" the complex poles, and the 


line 8p=6z which is the "second order" line. The important factors 


of these divided zones are as follows: 


(1) 


(2) 


(3) 


ep<90° implies a negative slope to the line which is used 
for computation as described in Section 5. 6p290° implies 

a positive slope, 

@prez implies that Sin(@z-@p)<O, and Sin(@z-@p) is the 
y-axis intercept of the line described above. O6p<@z implies 
that Sin(@z-0p)>O. If Sin(@z-@p)¢Othen the negative half 
cycles are plotted on the computational diagram, i.e. the 
values of the Sint for T=7A27, 37 4041, S0Altetc. 

If Sin(@z-@p)70 then the positive half cycles are plotted 
ise. the values of Sin? for T2=0% 1, 2A ZTE. (See 
Section 5). 

If the slope of the computational line Sin(@z-@p) <a macy 
is negative then the line must intersect all the half cycles 


of Sin T and an infinite number of extremal points exist, If 


the slope is positive then it intersects the function SinT 
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Fig.ll: Regions of Interest on Op vs. @z Plot. 
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(4) 


(5) 


at only a finite number of points and a finite number of 
extremals exist (this finite number may be zero). 

Op 7ez implies that the value T+@z-Op lies to the left of 
the value T on the computational diagram, If this 
condition is coupled with the condition that the slope of 
the computational line is positive then the only manner 
in which the equality expressed by Equation II can exist 
is for the computational line to "pierce" the curve of 
Sin T (the Sin T curve is concave downward; therefore, 
"piercing" means to intersect on the positive slope side and 
go "inside" the curve.) Any line which "pierces" must 
"emerge"; therefore, the extremals always occur in pairs. 
If the first extremal is called an "overshoot", then it 
can be stated that every "overshoot" is accompanied by 

an “undershoot"; however, it is possible that the time of 
the "overshoot" and time of the "undershoot" may coincide, 
@p <@z implies that the value T+@z-@p lies to the right 
of the value T on the computational diagram. If this 
condition is coupled with the condition that the slope 

of the computational line is positive then the equality 
expressed by Equation III can exist for a "piercing" 
condition, or it can exist for the condition in which the 
computational line approaches just to the horizontal 
distance Op-0z from the curve Sin T and no closer. If the 


latter condition occurs then an "overshoot" and an "undershoot' 
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f occur at the same time, defined by the particular value 

of T, and no extremals will occur after that t ime;the 

Slope of the time response goes to zero, but does not 

change sign. The first extremal always exists, but 

extremals after that time occur in pairs. 

These factors are of the utmost importance, and if they are 

not thoroughly understood the reader should review the procedures 
outlined in Section 5 for computing Tp. These factors are graphically 
illustrated in Figures 12, and 13. Figure 12 is a sample computational 
diagram for @p<0z. The values of Sin Tfor positive half cycles are 
drawn. One sample computational line has been constructed for the 
case @p<90° (@p=73° for the example). The points where Equation II 
is satisfied are depicted by short vertical lines on the computational 
line, Four sample computational lines have been drawn for the 
case @p790° (@p=108° for the example). Each line has been drawn for 
the same ©p but a different value of @z-Op. For line "A" the 
extremals exist as overshoot-undershoot-overshoot and no further 
extremals. For line "B" the first overshoot occurs, but the times 
of the first undershoot and the second overshoot coincide, and no 
further extremals exist. For line "C" only the first overshoot exists 
with no further extremals. For line "D" the value of 6z-0p has reached 
its maximum of 72° (@z2 maximum is 180° and @p is 108°), and then 
the first overshoot occurs at TPO; no further extremals exist. 
Figure 13 is a sample computational diagram for @p76z. The values 


of SinTfor negative half cycles are drawn. One sample computational 
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line has been constructed for the case @p<90° (@p=73° for the 
example). Four sample computational lines have been drawn for 
the case @p>90° (@p=108° for the example). For line "A" two 
overshoots and two undershoots exist. For line "B" one overshoot 
and undershoot exist. For line "C" the times of overshoot and 
undershoot coincide. For line "D" no extremals exist. The 
properties of the regions on Figure 13 may now be stated in terms 
of these factors: 
(A) Region I, Referencing (1), (2), (3), Line E-Figure 
12: All extremals exist, slope of the computational 
line is negative, and positive half cycles are 
plotted on the computational diagram, 
(B) Region II, Referencing (1), (2), (3), Line E-Figure 
132 All extremals exist, slope of the computational 
line is negative, and negative half cycles are 
plotted on the computational diagram, 
(C) Region III, Referencing (1), (2), (3), (4), Lines 
A,B,C,D =- Figure 13: A finite even number of 
extremals exist, slope of the computational line 
is positive, negative half cycles are plotted on 
the computational diagram, In region IIIA, the 
computational line will always appear as line D, 


i.e., there can be no overshoot in this region. 
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Diagram(6p<@z). 
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(D) Region IV, Referencing (1), (2), (3), (4), Lines 
A,B,C,D - Figure 12: A finite odd number of 
extremals exist, slope of the computational line 
is positive, positive half cycles are plotted 
on the computational diagram. In region IVA, 
the computational line will always appear as 
line C or D, i.ee., there can be only one over- 
shoot and no further extremals in this region. 
Regions IIIB and IVB are of special interest because they 
bound the line @z=@6p. Along this line the pole and zero of the 
system cancel each other and the system acts as a pure second 
order system. An unsuccessful attempt at cancellation would 
cause the system to operate in regions IIIB or IVB. Consider 
first the region IIIB and refer to Figure 13 and the computational 
lines for 6p 790°, Assume that @p=108° and @z is to be varied from 
@z=108° to @z=0° which will be a movement along the horizontal line 
@p=108° through regions IIIB and IIIA. Obviously when @z=@p, 
Sin (@z-@p)=O and all extremals exist and occur at multiples of 77 
(this is the second order system). As @z decreases the computational 
line translates upward rapidly; when 9z=106,3° [Sin(@z-0p)= -.03 | 
the point has been reached where the second overshoot and second 
undershoot occur at the same time (Tp=10.6rad). When 92=96.5° 
[Sin (0z-Op)= areal the point has been reached where the first 
overshoot and first undershoot occur at the same time (Tp=).5rad)3; 


this point is the dividing line between region IIIB and region IIIA. 
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For @2€96.5° no overshoot can exist. Consider next the region 
IVB and refer to Figure 12 and the computational lines for 
8p?90°, Again assume @p=108" and @Z is to be varied from 
@z=108° to @z=180° which will be a movement along the horizontal 
line @p=108° through regions IVB and IVA, As 92 icagieaiyr the 
computational line translates upward; when @z=113.7° | Sin (@z-8p)= 
0.1] the point has been reached where the first undershoot and 
second overshoot occur at the same time(Tp =7.35 rad); this point 
is the dividing line between regions IVB and IVA. For 627113.7° 
only one overshoot can occur. When @z=180° [ Sin (@z-6p)= 95 | 
The bY eeshoet occurs at its minimum time of Tp=1.57. It is 
possible to continue the reasoning in this fashion for various 
values of @p and thus summarize the properties of regions IIIB 
and IVB. This summary is presented in Table II. 

The next factor of interest for the third order system with 
One zero is the normalized magnitude of the extremals. This 
information cannot be plotted on a two-dimensional plot since 
Mpt is a function of Op, @z, and v. It is not difficult however 
to predict the order of magnitude of the overshoot. Along the 
line @p=6z the system is second order and the normalized magnitude 
of the overshoot can be read from the second order curve, Along 
the axis @p=0° the system operates as second order with one zero 
and the normalized magnitude of the overshoot can be read from 
Figure 6, Along the axis 6z=0° the system operates as pure third 


order and the normalized magnitude of overshoot can be read from 
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Figure 8(a). It is also possible to plot the line along 
which Mpt is unity. This is accomplished by returning to 
Eq'n IIIa and equating to final values 
Mof=f= 1+ SintBe1w) 9 WHT | sual Sin (tpp99-Gp)- Sin ("2 +0 2-9 +7x)] 
fe oe S'9 (19+ 02-02) = Sin (19+ 02-997) 
Sinope Sia (Tp 402-0) = Sin (0917) S19 (Tp + O2- GIT) 


Cos (Te 702) = Cos (PIPa tery) 
Which rs ods fPOOT When i 
7? = —-(024m7) 22°97 
This curve can be plotted byselecting a particular Y and plotting 


points on the 7p curves by the equation: 


G2= 297 ~0— Fe 
Curves have been plotted on Figures 10a and 10c, for selected 


values of 7, depicting the curve along which Mpt=1 (first 
overshoot in Figure 10a and second overshoot in Figure 10c). 

Note that these curves end at @z=p=180° - 7 which represents 

the origin of the S-plane. Beyond this point the pole and zero 
of the system fall in the right half plane, and the system is 
unstable. A limiting situation in determining magnitude of the 
extremals is encountered by letting 9z2=180-% (zero at the origin). 
When this occurs Sin (@z +% )=Sin 180° =O and the expression for 
Mpt becomes infinite, but the final value 5 E38 goes to zero so 
that the expression is meaningless. To evaluate the overshoot 


the value Z is brought inside the parentheses in Equation III: 
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f Vener 


[2272 7/7 (ue/for 82-82) = Sin (wd 40 + -epitls 


a 
on 2S (9277) 7 
X dp) a jz Ps 


SAT 5$:0(9217T) 


7 ; $3 P 
r Chie SS ’ eo JA 40 -B as 
ln? Zz 3nT Sindy e EAPe sm lomtlp 102 %) Sin (uplle 163-07 }é 


For ie Zero ve/, Wie Origin 


Z20 5 O2 = 180-T 5 S170 O2= SAT: 


Vor 
we X (4p) = Ge = dof te "[snee Sin (whe 1'80-4-07 ) = Sin (wate +180 -sp)]¢ 


Sin (opt7r) 


- = j Soleo tr 0 A ny | £022 Sr (wal4j0 “Op -r) Syl (we tye -o/\§ 


“pn 51n(optt) 


This equation is good along the line @z=180°-F until the point 
@p=180°- F is reached. Bringing P inside the parentheses? 
= - uh yr Sia 
Xp =e Yor S SinBP Sin(ugtp-Op-4)- S20 ID 5, tyy]e-Op) 
— P Sin? ¥ 


Psa? 
wolf " % 
-f_ ont Si02G2 Sin (wolip-Sp-@) _ ens Sin (uote -S)§ 


Wn® Snr Sia (6e17T) 
Buf e Oz= 180-7 , E ‘A, wz pede bx ra 
Sin (\80°-7- Ope (Phen 6p = Sin (udp 4180-7 -o7) 
ewe feebp =» Sin (wolle-p-® 


Sin t(epdtr 


Si2Pf s,_ (wasp ap) 
SAt Or 


Per : 
B Xg)2 eA ‘|- sin BBY owt hanse 


Wy4 Ae 


he Hrag op = )80°° TF 
= =k gunner le e walt /pon of — = hyd 2 = fe) 


ly,* Sin 


= | 2 +e : ad (Ginko yes 7) 


which is the expression obtained for a pure second order system. 
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Another limiting expression occurs when @p=180- 7 


(pole at the origin). When this occurs Sin (@p+7 )=Sin180°=0, 


but the final value i£55 becomes infinite so that the expression 


is meaningless. To evaluate Mpt the value P is brought inside 


the parentheses as before: 
watt ne 


X (dp) = EE, io » Sin (opie) e |s2 Si n(wya rp brea ep) +5:0laupr97-ap 2) § 


Psat Sn (0210) Sin (9d 1) 


‘a 517 OP 


Wan Sin (OPt r) 


lee 
pj) = rm ae Ta ee) [em Sn (tty Jo +82 -Op) t Sin (4tr40 109 -Sp i] § 
= alld his oe an To y 
op 4Phen - wl ber, by2 
ato So BP 7 SIN Of? 
Ly fuser» fae orden? [ene aE Oe PE + Sin (at Og apis 


a Gp = IZO- «as SinPa= ane 5 fenOpe-fenT t 


Xt) Oe eas bho Sin (wide 402420) 
oa re lectern 
(p= “Wps © Sin (02 +7) 
This equation is good along the line 6p=180° - 7 until the point 
@z=180° -7 is reached. Bringing Z inside the parentheses and 


oso? ter) 
Sar ena 


X tp) = Ee. oon Sia (wad 102427) 
; S10 Oe 


hedbng O22 /80-F s sin bpeSnT: 
« weld 
X (4p)= in e WA fen TF 
Bud wllp= 7 : 
arn no 
a Xcite) = Zen e He 


It may be noted that this is the expression for magnitude of the 


recognizing thats += 


first overshoot of a second order system with the exception of the 


additive constant im 
The character of the variation of the magnitude of the first 


peak may be summarized by means of a diagram as depicted in Figurel). 
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The coordinates are the @z and 6p axes used for the Pp plots, 

The line @z=@p has been drawn in. A value of 7 =60° (cos =.5) 

has been selected arbitrarily as an example. The value of 

Mpt=1.16 has been determined from the second order equations for 
placement on the @z=@p line. From Figure 6 values of Mpt for various 
@z have been placed along the @z axis. From Figure 8(a) values of 
Mpt for various @p have been placed along the 6p axis. The curve 
for Mpt=1.0 has been copied from Figure 10a (all values of Mpt 
assuming final value =1). The curves outlining the regions of no 
overshoot and of one overshoot only have also been drawn in. The 
lines @z=1809°-f and @p=180°% 7 have been drawn in to depict the 

right half plane (RHP) and left half plane (LHP). The values of 
overshoot along these latter lines are determined by the equations 
derived above, 

Approximate curves for Mpt=1.1 and Mpt=1.3 have been sketched in 

to show the form of the curves, This form is that the curves are 
convex toward the @z=@p line and all curves intersect at the origin. 
The value of Mpt becomes unbounded for values close to the 6z=180-7 
line and approaches zero along the @p=180-7 line. Mpt is positive 
in regions I and IV and negative in regions II and III. Itisa 
Simple matter to sketch a diagram of this sort for any value of T 
and thereby predict the action of the system both in time and 
magnitude, Similar diagrams for the first undershoot and second 


overshoot may also be quickly sketched, 
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No curves have been plotted for rise time (tc) since the 
computation is straightforward and the information already 
plotted gives a complete picture of the response of the 


system. 
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10, A New Look at Some Old Concepts: Dominance and Cancellation. 

It has been stated repeatedly in this thesis that the most 
common method of treating the higher order systems is by making 
a second order approximation. This is sometimes done by forcing 
the real singularities far out on the negative real axis (dominance) 
or by cancelling with a pole and a zero (cancellation). Each of 
these methods is performed in a rather vague manner; i.e., how far 
is "far out" in the LHP in the case of dominance or how close is 
"close together" in the case of cancellation. Another question is: 
what parameter of the response is of the most importance in making 
the second order approximation? Is the designer/analyst interested 
in time of response, magnitude of fluctuations, time of extremals, 
etc.? It is possible, by making use of the methods of this thesis, 
to take a new look at these concepts and perhaps make them a little 
less vague, 

Consider the case of a third order system with a troublesome 
pole on the real axis. Using the concept of dominance a rule of 
thumb says to make the pole ten times as great as the real part of 
the conjugate roots, but an examination of Figures 7 through 9 
discloses that making @p=30° causes a variation of no more than 
10% in the magnitudes of the extremals (the magnitudes are less), 
no more than 15% in the times of the extremals, and no more than 
30% in the rise time (the times are greater). While some or all 
of these variations may be too large, two points are demonstrated: 


(1) If the time parameters of the response are the important 
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factors then Op must be made smaller than if the 
magnitude parameters are the important factors. 

(2) The required distance of the pole from the origin varies 
with the value of 0%. For example, in the above case, 
if 7=30° then it would only be necessary to make the 
pole twice as far out as the real part of the conjugate 
roots; however if % =60° the pole must be four times as 
far out as the real part of the conjugate roots. 

It is fairly certain that the rule of thumb value of "ten times 
as great" will make the pole insignificant, but it is equally 
certain that it is not necessary to be so drastic in all cases, 
Utilizing the methods described here and armed with the require-= 
ments of the situation at hand, the designer can develop the most 
efficient configuration, 

Consider now the same system as above, but assume that the 
troublesome pole is to be removed by cancellation with a zero, 
Examination of Figures 10 and 1h show immediately that the require- 
ment of how close the zero must be to the pole to effect 
cancellation is dependent upon the location of the pole with 
respect to the conjugate roots, value of the angle 3, and the 
relative importance of time and magnitude. For @p<70° 6z can be 
as much as 15° either side of the value of Op and the time of 
extremals will change no more than 10%, and if wv < 70° the change 
in magnitudes will be on the order of 10%. On the other hand, 


ce @p>90° the value of @z becomes critical. This is best illustrated 
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by example: 


Given: 7 =60°; @p=110° (See Figure 16) 


8z Tp Mpt 

2.4159 2.8 1.958 

112° p50 e853 

er 3.14 1.160 (exact cancellation) 
108° Bee 1wAOB7, 

105° 3.6 908 

100° 4e55 819 

100° (no extremal ) 


For a value of P=,.3 this range of @z represents a variation of Z 
from .199 to .489. It is readily apparent that a very slight error 
in cancellation can cause a serious deviation from second order 
operation, It may also be noted an error of 9z>@p causes a 

smaller change in Tp and a greater change in Mpt than an error 

of @z<@p, Again, it becomes easy to determine the margin of error 
allowed to effect the desired response. A quick glance at the 
charts will determine for the designer his tolerance and which 


direction to weight the uncertainty. 
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ll. Conclusion... 
In order to simplify the analysis of the third order equation 
with one zero the time response equation has been developed in 


terms of the angles at the singularities: 


» wd, brn utllerr 
Xa? Ae ‘4 + singesin (ed o We?” saga g Soult 6a p10 


Ine Sin (G27) Say sin (6249) 


This equation has been manipulated to provide expressions for three 
time response parameters. 


(1) Times of extremals ( tee . 
po / fen ep 


~ Wel 
Sin (02-Op)€ = Sin (ath pn + 92-9) 
(2) Magnitudes of extremals (xX Ud) 


Xp) = EE, $, Lt 519 (6247) oX wotdye) Jon * eee 7? (wade 169-Pp)eSin (udp wail 


Sin S10 (0840 
(3) Times when response crosses the final value ( 4, ): 


2 ~ WdtEen o 
S ya Sin (62-Op)e [Es Sin liijlen +02 go +7) 
A simple graphical method has been developed and demonstrated for 


solving the transcendental equations above, 

The equations above have been utilized to analyze four types 
of systems: 

(1) Pure second order 

(2) Second order with one were 

(3) Pure third order 

(4) Third order with one zero 
Examples have been worked and a variety of curves developed to 
enable the designer and/or analyst to make a rapid appraisal of 
the response of the various systems, 


An analysis has been presented on new methods for achieving 


dominance or cancellation and evaluating the effects of errors in 


such attempts, 
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It is felt that a new insight has been gained into the 
actions of the third order system with one zero. With the 
methods and charts of this thesis it is possible to predict 
immediately the form of the response, and only a short time is 
required to compute the parameters and sketch the response accurately. 
For the designer interested in improving the reSponse of a second 
order system by adding real singularities, the means is now 
available for doing so exactly without resorting to the trial 
and error methods. 
There is a possibility that the methods of this thesis might 
be extended to other studies. Some that have occured to the 
author are listed here in the hope that others may be persuaded 
to attempt to see the methods: 
(1) A generalized method for expressing the time response 
of higher order systems in terms of the angles at the 
Singularities, 

(2) Combination of the method with that of Mitrovic to 
develop generalized curves for the time response on 
the Mitrovic Plane. 

(3) Extension of the equations of this thesis to develop 

charts for frequency response parameters. 

(4) Combination with methods such as those of Ross-Warren 

and Pollack to develop improvements and/or simplifications 


in the methods of compensation calculations. 
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